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1. Introduction 



Understanding the quark and gluon structure of the nucleon is a vital component of our en- 
deavor to understand how QCD gives rise to the properties of the observed universe and is the focus 
of frontier experiments in nuclear and particle physics. Lattice QCD provides a unique tool to study 
nucleon structure from first principles, and as summarized in the plenary taUc at this conference by 
J. Zanotti [|l]] many successful techniques have been developed to calculate form factors and gen- 
eralized form factors of the nucleon. With the development and availability of Petascale computer 
resources, we are now entering an era in which high precision calculations of nucleon structure are 
becoming feasible, which necessitates a fresh examination of statistical and systematic uncertain- 
ties. Here we describe the results of several developments in our collaboration's efforts to enter a 
new regime of precision in calculating hadron structure. The results in this talk will focus on calcu- 
lations with a hybrid action combining domain wall valence fermions with improved staggered sea 
quarks, as described in Ref. ||^, and that utilize the extensive set of configurations with dynamical 
improved staggered quarks generated by the MILC collaboration [^. We will also refer to recent 
calculations using dynamical domain wall configurations [Q, |5|] described in the talk by S. Syritsyn 

We emphasize some technical aspects of the nucleon form factor calculation which are neces- 
sary for good controls of both statistical errors and lattice artifacts. A more thorough description 
of the calculation details and simulation parameters of our mixed action project can be found in 
Refs. [Q] and [^. This proceedings is organized as follows: In Section ^ we present an analysis 
which estimates the increase of the numerical cost as we go to the physical pion mass. In Section || 
we discuss an improved technique in the calculation of the backward propagators for the nucleon 
three-point correlation functions. Studies of correlations in the lattice data are given in Section 0. 
We discuss possible finite volume effects and discretization errors of our calculations in Section ^. 
We show selected form factor results and comment on the QCD evolution of our calculations of the 
flavor non-singlet nucleon angular momentum in Section ^ followed by conclusions in Section^ 



2. Numerical Cost of Precision Calculations 



As we progress to ever lighter pion masses, it is important to quantitatively estimate the statis- 
tics required to achieve a specified level of precision as a function of pion mass. From the per- 
spective of fluctuations produced by a set of gauge configurations, the three-point function and 
corresponding two-point function with the same source sink separation behave similarly, so the 
two-point correlation functions C2{t) = {J{t)J(0)), where J denotes the nucleon source, are a use- 
ful measure of the statistical fluctuations. Since the variance in C2{t) is generated by a source 
containing three quarks and three antiquarks, it receives contributions from both three-pion and 
two-nucleon states, and the corresponding signal to noise ratio is therefore given by 



Signal (-/(O-^W) 



Noise -^V(l^(0^(0)P)-((/(0/(0)))2 -^VBe-W-Ce-^^~' 
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Figure 1: The left panel shows the signal-to-noise ratio of the nucleon two-point correlation function on 
Asqtad lattices at nij^ ^350 MeV, where beyond t = 12, the exponential decay is given by {Mn — \niji). The 
right panel shows the corresponding "effective mass" of the signal-to-noise ratio in (a), where the horizontal 
line is the measured — jm-n of the ensemble. 




Figure 2: The left panel shows the extrapolation of the signal-to-noise ratio to the physical pion mass. The 
right panel shows the exponential increase of the number of configurations needed to maintain 3% accuracy 
at separation t= 10 as one approaches the chiral limit. 



which decreases exponentially with {Mn — ^mjt)t. Figure |l(a)| shows a typical result for the signal 
to noise ratio as a function of the time separation, which displays the expected exponential decay 
in Mn — ^niji at large t. This is clear from Figure 1(b) , where we show the "effective mass" of 
the signal to noise, compared with the measured Mn — ^niji of the ensemble, which is denoted 
by the horizontal line. Using comparable calculations for the three lowest pion masses, ~ 300, 
350 and 500 MeV, we can extrapolate the signal-to-noise ratio to the physical point, as shown 
in Figure 2(a). Correspondingly, the number of configurations required to attain 3% accuracy is 



shown in Figure 2(b) 



3. Coherent Sink Techniques to Increase Statistics 

Given the need for 5,000 to 10,000 independent measurements to overcome the exponentially 
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m 



,DWF 



[MeV] # confs # meas 



20^ X 64 0.007/0.05 
20^ X 64 0.01/0.05 
28^ X 64 0.01/0.05 



293 
356 
353 



464 3712 
628 5024 
274 2192 



Table 1: Numbers of measurements for the three lowest pion masses. 



increasing noise produced by three-pion states as the pion mass is decreased, we generate 8 inde- 
pendent measurements on a lattice of time extent T = 64 as follows. On the first configuration, 
we place sources at space-time positions (0,0), (L/2, 16), (0,32), and (L/2,48), and calculate 12 
sets of propagators which we will refer to as forward propagators in the usual way. Using the for- 
ward propagators from the source (x,, 7]), we create a momentum projected nucleon sink at time 
To away at {xi, Tt + To). A conventional calculation would combine forward propagators from the 
source at {xi, T,) and backward propagators from sink at (x,-, 7]- + Tq) to obtain the relevant two-point 
function, requiring 4 sets of inversions to treat all 4 sources. In contrast, we calculate a single set 
of coherent backward propagators in the simultaneous presence of all 4 sources. Combining these 
coherent backward propagators with the forward propagators from the f'^ source yields the physical 
result for the i'^ source with the i''' sink plus terms that vanish by gauge invariance when averaged 
over an ensemble of configurations. In addition, on the same lattice, using the forward propagators 
from the source (x,, 7]), we also create a momentum projected antinucleon sink a time Tq away 
at (x,, 7] — Tq), and perform an analogous calculation for coherent antinucleon propagators. It is 
straightforward to relate the matrix elements of our twist-two quark operators in an antinucleon to 
the desired results in a nucleon. The net result is that, given a set of forward propagators, we obtain 
eight measurements at the cost of two rather than eight sets of inversions. 

To minimize correlations, the sources on the next configuration in the ensemble to be analyzed 
are located at space-time positions (L/2,0), (0, 16), (L/2, 32), and (0,48), and subsequent config- 
urations are each shifted by a displacement L/2. The independence of these lattice measurements 
is addressed in the next section. We used these coherent sink techniques on the three lowest pion 
masses, the parameters of which are given in Table ^. 

4. Statistical Analysis 

4.1 Binning and Autocorrelations 

A crucial question concerning our calculations with 8 measurements per lattice is the statis- 
tical independence of measurements within a single lattice and between subsequent lattices. One 
standard test of correlations is binning potentially correlated measurements and observing the de- 
pendence of the jackknife errors on the bin size. Figure |3| shows the results of measurements with 
five different binnings for the two point function C2 measured midway between source and sink 
(t = 5) and at the source-sink separation (t = 9), and three current operators, Jx, Jy, and 7,, mea- 
sured midway between the source and sink (t = 5) measuring the electric form factor, Ge, the 
magnetic form factor, Gm, and Ge, respectively. Bin size 1 treats each measurement separately, 
size 2 combines nucleon and antinucleon from the same source, size 4 combines two nucleon and 
two antinucleons, size 8 include all 8 nucleons and antinucleons on a single lattice, and size 16 
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C2(5) 



lmJx(5) 
~GE 



Re Jy(5) 
~GM 



Re Jt(5) 
~GE 



Figure 3: Jackknife errors with bin sizes ranging from 1 to 16 measurements as described in the text. 
Independence of the error with bin size indicates negUgible correlations between measurements. 
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(b) Tint for three-point functions. 



Figure 4: Representative integrated autocoiTelation time for two-point and three-point correlation functions. 
For completely decorrelated measurements, the integrated autocorrelation time is normalized to 1/2. 



combines two sequential lattices. As is clear from Figure ^ the negligible change in the errors with 
bin size indicates negligible correlations. 

One can also check the correlations between the measurements by calculating the integrated 
autocorrelation time [H], defined as 



Ncu, 

Tint = 1/2+ £pW/p(o), 

n=l 



(4.1) 



where p{n) is the autocorrelation function between «th and 0th measurements. The separation 
between two independent measurements would be iTjnt. 

We calculated Tint in our measurements by treating each source as an individual measurement. 



Figure 4(a) shows the result for a zero-momentum projected two-point correlation function at a time 
separation t = 5. The horizontal axis is the cut in the summation for the integrated autocorrelation 
time in terms of the number of measurements, which is 8 per lattice. One can see that Tint reaches a 
plateau at around 0.7, meaning the measurements are already de-correlated for every other sources. 
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(GeV^) 
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(a) 



(b) 



Figure 5: Example of overdetermined analysis for the F2 form factor as a function of including large 
boosts of source and sink. The left panel includes the source momentum components (-2,0,0) or larger, and 
sink momentum (-1,0,0). The right panel has had these data points removed. 



A similar result for the three-point correlation function is shown in Figure [4(b) , in which Tjnt is 
very close to 1/2, indicating that the measurements from adjacent sources are independent. These 
results are consistent with the binning study as discussed previously. 



4.2 Correlations Among Different Momentum Transfer 

To maximize the hadron structure information determined from a given set of lattice config- 
urations, we use the overdetermined analysis introduced in Ref. to simultaneously extract a 
specified set of generalized form factors from as many different combinations of twist-two opera- 
tors and source-sink momenta as possible. One potential liability of this approach is the admixture 
of noisy measurements arising from source and sink momenta that are sufficiently high that the 
data are subject to large statistical errors. If such data are obviously consistent with the accurate 
data, they do not affect chi-squared fits, and there is no significant bias in removing them from the 
analysis after the fact. However, frequently the offending data gives the superficial appearance of 
being statistically inconsistent with the accurate data, if one ignores correlations. A typical exam- 
ple is the measurement of Fj" on a 28^ x 64 lattice at niji ~ 350 MeV, as shown in Figure 5(a) , 
where it appears that 6 out of 21 data points lie significantly above the reference curve. 

The question is whether 6 independent measurements lie several a off the fit to the accurate 
data, which would be a highly statistically significant discrepancy, or if the data are highly corre- 
lated so that only one or two degrees of freedom have fluctuated randomly. Hence, we calculated 
the correlation matrix 

(4.2) 



'"'7 ~ ^i}/ \l^n^ib 



where the covariance matrix is defined as 



(4.3) 



M=0 



p'f^ is the jackknife sample of the Z^" momentum (hence the factor — 1), and F, is the ensemble 



Jh : 



;th 
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average of the momentum transfer. The resulting correlation matrix for the six outliers is 



(4.4) 



/ 1 0.822855 0.654945 0.639426 0.552542 0.520037 \ 

0.822855 1 0.874902 0.649138 0.712735 0.742609 

0.654945 0.874902 1 0.562507 0.697565 0.701146 

0.639426 0.649138 0.562507 1 0.818721 0.408017 

0.552542 0.712735 0.697565 0.818721 1 0.584403 
V 0.520037 0.742609 0.701146 0.408017 0.584403 1 / 

which shows very strong correlations between the data. 

To quantify the correlation, we calculated the of the outliers relative to a reference curve, 
which involves two steps. First, we defined the reference curve by performing a correlated 
dipole fit to all the data points (except for = 0), including the six "outliers". Second, the x^ of 
offending data points was calculated by 



n=0 



(«) 



(4.5) 



where Fi is the expected form factor result at momentum from the fit, and A'^ is the number of 
outliers. If x^ is close to 1, then the data points are consistent with the fit curve, and hence, the rest 
of the data points. In this example, the x"^ is found to be 1 .9 ib 1 . 1, indicating that the deviation 
is not statistically significant. In contrast, a naive visual analysis of the curve ignoring correlations 
would lead to the erroneous conclusion that uncorrelated points are away from the curve by two a, 
in which case the X^ 1^ would be 4. This justifies the systematic exclusion of large source and sink 



momenta that give rise to noisy outliers, the result of which is shown in Figure 5(b) 



5. Control of Lattice Artifacts 

5.1 Finite Volume Effects 

As we go to lighter pion masses while holding the lattice volume fixed, the effects of finite 
volume become more and more important. While it is widely accepted that nijiL > 4 is necessary 
to avoid sizable finite volume effects, there are still controversies on how big the volume should be 
to have a good control over finite volume effects for nucleon physics. Here we do not attempt to 
address this question from a theoretical point of view. Rather, we present numerical evidence from 
our mixed-action calculations with two different volumes at a pion mass of roughly 350 MeV, to 
estimate how large the finite volume effects might be with this particular action and chosen lattice 
parameters. 

In Figure 5(a)| , we compare the results for the isovector Dirac form factor from the 20^ x 64 



and 28^ x 64 ensembles, corresponding to physical volumes of roughly (2.5fm)^ and (3.5fm)^, 
respectively. The solid curves are dipole fits to the lattice data with < 0.4 GeV^. The fit pa- 
rameters agree within errors, showing that there are no statistically significant differences between 
these two volumes. Since the slope of the form factor at 2^ = gives the Dirac radius of the nu- 
cleon, one can infer that the results for the Dirac radius obtained from these two volumes do not 



show significant finite volume effects. We present the comparison of the Dirac radii in Figure 6(b) 
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o = 356 MeV [L=2.5 fm] 
□ m = 353 MeV [L=3.5 fm] 
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(a) (b) 
Figure 6: (a) Isovector Dirac form factor from two different volumes, (b) Isovector Dirac radius. 





20^ X 64 (L w 2.5 fm) 


28^ X 64 (L 3.5 fm) 


Ay 




0.274(10) 


0.290(10) 


-0.016(14) 


gA 


1.161(14) 


1.153(14) 


0.008(20) 



Table 2: Comparison of results for (rj") and gA from two different lattice volumes. Ay is the difference 
between the L = 2.5 fm result and the L = 3.5 fm result. The error on Ay is calculated by adding errors from 
the two volumes in quadrature. 



where we also include the result at the pion mass of 293 MeV. And the star is the phenomenological 



value as obtained in ||10|]. 

To estimate the finite volume effects quantitatively, we calculate the differences, Ay, from the 
small- volume (2.5 fm) and large- volume (3.5 fm) calculations for the isovector Dirac radius, {r\^), 
and nucleon axial charge, gA- Since these two ensembles are statistically independent, we simply 
calculate the errors on Ay by adding errors in quadrature. The results are given in Table H We 
can see that the differences are statistically consistent with zero, suggesting that the finite volume 
effects are negligible compared to statistical errors. 

5.2 Comparison of Domain- Wall and Mixed Actions 

Since domain wall fermion formulation is automatically 0{a) improved, we expect the dis- 
cretization error in the full domain wall calculation to be small. In our full domain wall calculations 
of the nucleon form factors, this is found to be the case, as described in Ref. [^]. To assess the dis- 
cretization error of our mixed-action calculations, here we compare some of the nucleon structure 
results from both actions as available at the time of the lattice conference. For latest domain wall 
results with improved statistics, see Ref. 

In Figure 7(a) we show the isovector Dirac form factor from the fine {a^^ 2.346 GeV) 
domain wall and the coarse {a^^ 1.588 GeV) mixed-action calculations, both with a pion mass 
of roughly 350 MeV. Over the whole range of available to us, both actions give statistically 
consistent results. We note that the lattice scale of the fine domain wall ensembles used here is 
a crude estimate obtained in Ref. ||^. Using the Sommer parameter, the lattice scale for the fine 
domain wall enembles was found to be = 2.42(4) GeV [pi]. If the latter were used, the Q 
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(a) Isovector Dirac form factor. 
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(b) Isovector Dirac radius. 
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0.5 



(c) Nucleon axial charge. 



Figure 7: Comparison of some physical results from the mixed-action and domain-wall calculations. In (b) 
and (c), the a^' = 2.346 GeV domain wall results are slightly shifted to the right for clarity. 



values would shift to the right, and the pion mass would become larger. The combined effect would 
shift the domain wall results upward to be even more consistent with the mixed-action data, which 
is an indication that the discretization error for the mixed-action calculation is small. 



Figure 7(b) shows the isovector Dirac radius from the coarse mixed-action calculation and 
domain wall calculations at two different lattice spacings. The curve is a one parameter {B\q) fit 
to the three mixed action data points using the 0{e^) small-scale-expansion (SSE) chiral formula 
as given in [[l^], with the low energy constants fixed to phenomenological values. One can see 
that all the data points fall on the same curve, showing that the mixed action results qualitatively 
agree with the domain wall data. The same holds true for the nucleon axial charge, as shown in 
Figure 7(c) , where all the data show little pion mass dependence, and all lie on one horizontal line 
which is a few percent below the experimental result. 

Since in both the mixed-action and the fine domain wall calculations, the lightest two pion 
masses are comparable (~ 300 and 350 MeV), we are able to do a more quantatitive comparison. 
Table || gives the differences between the mixed-action (MA) and domain- wall (DW) results, A^, 
for {r\^) and gA, at these two pion masses. Once again, the differences are found to be consistent 
with zero within errors. 
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(rf ) [fm2] 


gA 




mj, ^ 300 MeV 


« 350 MeV 


300 MeV 


« 350 MeV 


MA 


0.311(22) 


0.274(10) 


1.176(32) 


1.161(14) 


DW 


0.321(14) 


0.270(17) 


1.150(30) 


1.162(30) 




- 0.010(26) 


0.004(20) 


0.026(44) 


-0.001(33) 



Table 3: Estimates of differences between mixed-action and domain wall calculations for {r\^) and gA- 
is the difference between the mixed-action and domain-wall action results. The error on A„ is calculated by 
adding the errors from the two actions in quadrature. 




1 

0.2 0.4 0.6 0.8 1 1.2 

[GeV^] 



Figure 8: Isovector electric form factor from the mixed-action lattice calculations along with the Kelly 
parametrization [|4p of the experimental data. 



6. Select Physical Results 

6.1 Electromagnetic Form Factors 

Complementary to the results for the nucleon electromagnectic form factors mentioned in 
previous discussions, in Figure |8| we show the isovector electric form factors from several masses, 
including results from previous calculations [ p^ ] at pion masses of 495 MeV, 597 MeV and 688 
MeV. To avoid duplication, we do not show results with the pion mass of about 350 MeV, which 
was already discussed in Section |5.1[ The curves through the lattice data are one-parameter dipole 



fits. The bottom solid curve is Kelly's parametrization of the experimental results []14[]. While 
in previous sections, we have seen that the chirally extrapolated lattice result for the Dirac radius 
is compatible with the experiment, here we see that over the range of momentum transfer up to 
< 1.2 GeV^, the lattice results show a monotonic decrease towards the experimental curve as 
the pion mass gets smaller. To have a direct comparison with the experiment for the dependence 
on the momentum transfer, we need to perform a chiral extrapolation with rriji and dependences 
taken into account simultaneously [p^, which is still a work in progress. 
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Figure 9: The left panel shows the quark spin and orbital angular momentum contributions to the nucleon 
spin from ref.[||]. The right panel show evolution of the quark angular momentum with respect to the scale 



6.2 Quark Angular Momentum Contribution to the Nucleon Spin 

The origin of the nucleon spin is a forefront question in contemporary experiment and theory. 
Figure ^ shows our recent lattice calculation of the quark spin and orbital contributions to the 
nucleon spin, which makes a strong case that the sign of the spin contribution of a given flavor, AT 
and the sign of the orbital contribution L for the same flavor are opposite at the scale 4 GeV^ at 
which the lattice calculation was renormalized. 

Since this behavior is strikingly different from simple familiar models, it is interesting to 
consider its origin and significance. In a mean field model, in which a quark satisfies the Dirac 
equation in a central potential, for a nucleon with spin projection in the positive z direction, the 
upper component is an ^-state with spin up and the lower component is a P-state with orbital 
angular momentum projection -i-l and spin down. Thus, generically, the spin and orbital angular 
momentum are aligned in mean field theory. Clearly, the reason we solve QCD from first principles 
is to go beyond models, and there is no reason that the nucleon should necessarily be consistent 



with mean field theory. However, stimulated by a recent paper by Thomas []15|], it is interesting to 
note that the sign of L"^'' is strongly scale dependent. 

We restrict our attention to the flavor non-singlet sector, for which there are no disconnected 
quark diagrams, and no mixing with gluons, and consider the evolution of orbital and spin contri- 
butions shown in the right panel of Figure ^. The spin contribution AL"^'^ is large and is conserved 
under QCD evolution as shown by the dashed line. However, at one-loop level the total angular 
momentum has the simple evolution given by its anomalous dimension 

where t = ln(-^). Because the largest contribution to the total angular momentum is the spin, 

which cannot evolve, the full change under evolution must arise from the orbital contribution L"^'^, 
which must therefore vary significantly with the scale. As shown by the solid line in the right 
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panel, the change required in L"~" is so substantial that it changes sign when the scale becomes 
low enough. Although surely one loop evolution is not quantitatively reliable and evolution below 
1 GeV is suspect, it is clear that even the relative sign between spin and orbital contributions is 
scale dependent, so there is no reason our lattice calculation at 4 GeV^ need to be consistent with 
simple quark models or mean field arguments. 

7. Conclusions 

We have investigated several aspects of the nucleon form factor calculations in an attempt 
to address some issues which are essential for achieving high precision calculation in the chiral 
regime. In particular, as the numerical cost increases substantially as our calculations move to 
Ughter pion masses, we employed the coherent sink technique in the calculation of the backward 
quark propagator to reduce the computational cost. To obtain reliable statistical errors, we also 
studied the autocorrelations in the measurements and found no strong correlations between mea- 
surements from adjacent sources. Our results at a pion mass of 350 MeV from two lattice volumes, 
(2.5fm)^ and (3.5fm)^, show no significant differences, suggesting negligible finite volume effects. 
We also compared the mixed-action results with the full domain wall calculations at a finer lattice 
spacing, and found that the discretization errors from using the mixed action is quite small. 
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